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We present a canonically invariant form for the generalized Langevin and Fokker-Planck equations. 
We discuss the role of constants of motion, and the construction of conservative stochastic processes. 
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Consider the usual Langevin equation: 



dV 

mqi = qi + £i(t) (1) 



with £j(t) Gaussian white noises < £i(t)£j(tf) >= 2TSijS(t—t'), and T the temperature of the thermic bath. Rewriting 
this as a set of phase-space equations: 



(2) 



we notice two things. First, the form of the first equation is restricted to a hamiltonian H of the form H = 
YliPi I 'm + V (<?) ■ Second, the interaction with the bath has introduced an assymmetry in the treatment of coordinates 
and momenta. We shall in what follows formulate the Langevin and Fokker-Planck processes in a way that treats all 
phase-space variables on an equal footing. 

In general, Langevin equations can be motivated jfl 0] by considering the system with Hamiltonian H, coupled to 
an infinite set of harmonic oscillators with random phases at some initial time and energies given by equipartition at 
temperature T. Upon solving for the oscillators, and reinjecting their dependence on the equation of motion, one gets 
a Langevin equation which can be made Markovian by a suitable choice of distribution of the oscillators' frequencies. 

Actually, d|,0) are associated with a particular coupling of the form: 



Hcoup ^ qi 



N 



(3) 



where y l a are the coordinates of the oscillators of frequencies uj l a . 

In order to obtain a canonically invariant generalization, one can repeat the exercise with a coupling with the 'bath' 
of oscillators of a more general form: 



H coup = J2J2^ A a G i(l,P)y i a + BiGl{q,p)yi] (4) 

i a—1 

Performing the usual steps, one arrives at the following Langevin equation, valid for any phase-space variable A(q,p), 
in particular the coordinates and momenta qi,Pi: 

A = k{A, H} + ]T{A G.H^t) + {Gj,H}) (5) 
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We have explicitated an inverse time-constant k. Here {A, B} — J2i (f^ff^ ~ are ^ ne Poisson brackets. 

The Gi(p, q) (i = 1, R) are R arbitrary phase-space functions, originally controlling the manner of the coupling 
to the bath as in Eq. The noise is white and gaussian as before. 

In fact, we do not need to go into the details of the derivation, because it will be shown in what follows that 
this equation is a bona-fide Langevin equation in that it leads to the canonical distribution at temperature T. The 
definition of Equation (|^) is completed by specifying that it should be understood in the Stratonovitch sense: in a 
discretized form all phase-space functions in the right hand side have to be evaluated as an average of their values in 
the previous and the incremented time. 

Indeed, one can adopt tto's convention so the r.h.s. is evaluated in the previous time, and the equation now reads: 

A = k{A, H} + J2{A Gj}(£ 3 (t) + {G 3 ,H}) + T{G„ {G 3 ,A}} (6) 



Two particular cases are Gi = — (ft V i ; re = 1 , which yields (|1|;|2|), and Gi — pi V i ; re = which yields the massless 
version of ([!]) . Note that in general the Poisson brackets between the Gi need not vanish, in which case the equations 
(|J) cannot be taken through a canonical transformation to the form (Q). 

Let us now turn to the (Fokker-Planck) equation satisfied by the probability distribution P(q,p,t). It is a simple 
exercise (see || or [Q) to obtain this directly from equation (||) or (|). The result is: 

BP 

— + re{P, H} = ^{G„ {G 3 ,H}P + T{G V P}} (7) 

3 

It is now clear that P — exp(— H/T) is a stationary solution of ([?]). Again, with the choice Gi = — V i ; re = 1 we 
obtain the Kramer's equation. If instead we make Gi = pi V i ; re = we obtain the usual Fokker-Planck equation 
for diffusion without inertia. 

By writing (A)(t) = J ' dq dpA(q,p)P(q,p,t) we obtain for the evolution of the average of an observable A(q,p): 



d{A)(t) 
dt 



({A, H}) ^2({G U A}{G h H}) - T^({G,, {A, GJ}> (8) 



All three equations (|^), (Q) and (^) are canonically invariant in form: a canonical transformation of variables is 
obtained directly by transforming H and the Gi. Furthermore, the pure Hamiltonian term and the bath-coupling 
terms are explicitly separated. 

Equilibration 

In order to study the equilibration properties, let us define an H-function as Q || : 

W(U) = J \q \ P V(q,p,U)(TlnV(q,p,U)+H(q,p)) (9) 

which cannot increase, since: 

< (10) 

i 

Note that only the bath-terms contribute. 

If the equilibrium measure exists, TL is bounded from below, and we have that 

{G t ,H}P + T{G u P}^0 V* (11) 

If we parametrize P(q,p, t) as: 

P(q,p,t) = Q(q,p,t)exp(-H/T) (12) 

the limit ( |ll|) implies that, once stationarity is achieved : 

{G u Q} = Vi (13) 

Using this equation, we have that Q = {H, Q}. Since at stationarity ( |l3| ) has to be valid at all times, we obtain the 
necessary conditions: 
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{G t ,Q} = ; {G t ,{H,Q}} = ; {d, {H, {H,Q}}} = ; ... (14) 

In the usual Langevin case (|[||), Gi = £j, H = ^2 i Pi/2m + V(x) and the first two sets of equations suffice to prove 
that Q — constant is the only stationary solution. 
Constants of Motion 

Suppose the Hamiltonian has some constants of motion {H, K a } = 0. Depending on the choice of Gi, these constants 
will be preserved or not by the coupling with the bath. Indeed, Eq. (|J) implies: 

^ =J2{K a ,G j m(t) + {G j ,H}) (15) 

3 

The evolution of K a is then purely dictated by the heat-bath, and will be 'slow' in the small noise limit. 

If we wish to construct a i4T -preserving noisy dynamics we have to choose the Gi such that {K a , Gi} = V i. 

If, on the other hand, we couple the system to the bath through some constants of motion, that is Gi = Ki, the 
Langevin dynamics for any A(q, p) becomes: 

A = K {A,H} + Y J {AK 3 mt) (16) 

i 

which expresses the fact that the system receives random kicks in the direction generated by the Kj . 

An extreme and rather amusing form of this is the case in which we put a single G = H and k = 0. We then have 

A = {A,H}at) (17) 

and the associated Fokker-Planck equation: 

^=T{H,{H,P}} (18) 

The system diffuses back an forth along its classical trajectories. The probability distribution tends for long times to 
the smallest invariant structure compatible with the original distribution, and the entropy J dqdpP(q, p, t) In P(q, p, t) 
becomes stationary. 

Motion within a Group 

Another simple application is the construction of a heat-bath dynamics on a group. Suppose the Hamiltonian is 
constructed in terms of the generators Li of a group, satisfying {Li, Lj} — CijiLi. Then, 

Li = {L h H} + CwLifo® + {Lj,H}) = Cy 7 (wj + + Ci. J iC jsr Lo s L r L l (19) 

where the 'angular velocities' are defined as u>i = dH(L) / dLi. The group invariant J^-Z/f is clearly a constant of 
motion. 



In summary, we have presented a manifestly canonical-invariant form of the langevin and Fokker-Planck equa- 
tions. Within this formulation several questions originating from the underlying classical mechanics become more 
transparent. 
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